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Abstract 

Using a stochastic cellular automaton model for urban traffic flow, we study and 
compare Macroscopic Fundamental Diagrams (MFDs) of arterial road networks 
governed by different types of adaptive traffic signal systems, under various 
boundary conditions. In particular, we simulate realistic signal systems that 
include signal linking and adaptive cycle times, and compare their performance 
against a highly adaptive system of self-organizing traffic signals which is de- 
signed to uniformly distribute the network density. We find that for networks 
with time- independent boundary conditions, well-defined stationary MFDs are 
observed, whose shape depends on the particular signal system used, and also 
on the level of heterogeneity in the system. We find that the spatial hetero- 
geneity of both density and flow provide important indicators of network per- 
formance. We also study networks with time-dependent boundary conditions, 
containing morning and afternoon peaks. In this case, intricate hysteresis loops 
are observed in the MFDs which are strongly correlated with the density het- 
erogeneity. Our results show that the MFD of the self-organizing traffic signals 
lies above the MFD for the realistic systems, suggesting that by adaptively ho- 
mogenizing the network density, overall better performance and higher capacity 
can be achieved. 

Keywords: macroscopic fundamental diagram, traffic signal system, 
simulation 
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1. Introduction 



A central goal of traffic science is the formulation of appropriate macro- 
scopic variables characterizing and relating demand and performance of road 
infrastructure. On the level of a single street ( or freeway), the fundamental 
diagram (FD), introduced by iGreenshieldsl ( 19351) . expresses flow as a function 
of density. Fundamental diagrams for a single link are generically unimodal, 
describing a free-flow regime at low densities and a congested regime at high 
densities^. It is far from clear, however, to what extent such simple relations 
should extend to more complex systems such as urban road networks. Early 
studies in this direction date back at least to Godfrey ( 19691) . The first con- 
vincing empirical evidence that congested urban networks can display simple 
relationsh ips between network- aggregated dema nd and performance was pre- 
sented in iGeroliminis and Daganzol (|2007l 120081) . These works clearly indicate 
the existence of a Macroscopic Fundamental Diagram (MFD) in the city of Yoko- 
hama, relating network-aggregated production and accumulation^. Analytical 
theories attempting to explain t he e xistence o f MF Ds have been developed by 
Daganzo and Geroliminisl (l2008h andlrlelbingl <|2009T). and match the Yokohama 



data quite well. The discussion in lDaganzo and Geroliminisl {2008) analyzes the 
effects on MFDs of applying different signal timings, by treating traffic signals 
as exogenous capacity constraints. One of the main aims of the current work is 
to study and compare MFDs in networks governed by different types of adaptive 
traffic signal systems. 

Given the existence of MFDs, it is natural to a sk under what conditions 
should they be observed? In previous work on MFDs, Daganzo and Geroliminisl 
(2008) postulate sufficient regularity conditions under which MFDs should be 
expected to exist, including slow-v arying and dis tributed demand, and ho- 
mogeneous network infrastructure. iHelbinel (|2009h argues that the details of 
MFD curves should be expected to depend not only on the aggregated den- 
sity, bu t also on the spat i al den sity distribution. Taking these observations 
further. [Mazloumian et aL ( 2010l) argue that the aggregated flow should in fact 
be a function of both the aggre gated density and also its spatial variation. 



Geroliminis and Sun TRBl (|201ll ) demonstrate using empirical data that while 



strict homogeneity of traffic states is not necessary to observe a well-defined 
MFD, the spatial distribution of density is indeed a key quantity. In this work 
we study the spatial heterogeneity of both density and flow, and demonstrate 
how their behavior can be used as predictors of network performance. 

In particular, we discuss the relationship between the time evolution of spa- 
tial heterogeneity an d hysteresis. Hysteresis has been observed and studied in 
freeway networ ks in Geroliminis and Sun TRAl ( 2011 ). While hysteresis was 
not observed in Geroliminis and Daganzol ( 20081) . a careful empirical study of 



1 Even in this simplest of cases, however, it is typically found exp erimentally th at significant 
scatter is observed in flow-density relations of congested links; see[kcriicr (199^). 

2 The production and accumulation are surrogates for the flow and density, and are more 
readily measured in empirical trials. 
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MFDs in the city of Toulouse was undertaken by iBuisson and Ladier 
in which hysteresis effects were clearl y observed. A theoretical explanation of 
th e clockwise hysteresis loops found in Buisson and Ladier ( 20091) was presented 



Gavah and Daganzol ( 201ll ). by treating the network as a dynamical system 
and performing a stability analysis. This analysis, which assumed a strict ho- 
mogeneity in the network, suggested that clockwise hysteresis loops should be 
typical, while anticlockwise hysteresis loops should be rare. If the constraint 
of perfectly uniform loading is relaxed however, different behaviour can arise. 
For example, for networks corresponding to commuter corridors, i.e. portions 
of an arterial network wedged between a strong source (e.g. a residential area) 
and a strong sink (e.g. the central business district), the input/output rates 
on the boundary may be significantly higher than in the bulk of the network. 
Other situations where such behavior might arise include arterial networks in 
the presence of perimeter control. We observe from our simulations that for 
such systems both clockwise and anticlockwise hysteresis loops can generically 
appear, and we show that these o bservations can be explained using a modified 
version of the model presented in iGavah and Daganz 

Our s imulations util i ze a s tochastic cellular automaton (CA) model, intro- 
duced bv lde Gier et al.l (|201lh . This model is mesoscopic, in the sense that al- 
though individual vehicles are modeled, fine-grained details of individual driver 
behavior are deliberately treated in a course-grained, statistical, manner. While 
details of vehicular motion through intersections are deliberately ignored, real- 
istic signal phasing at intersections is included in the model. In fact, the model 
was specifically designed to provide a simple and fast way to study arbitrary 
traffic signal systems, on arbitrary networks. Using this CA model we study 
the existence and shape of MFDs for three specific traffic signal systems, using 
both time-dependent and time-independent boundary conditions. In particu- 
lar, we simulate variants of the SCAT^I traffic signal system, which is cur- 
rently employed by numerous road authorities worldwide, including in Sydney 
and Melbourne. In order to study the effect of increased adaptivity on MFDs, 
we then compare these results for SCATS with the highly-adaptive (idealized) 
self- organizing traffic lights (SOT L) system, originally introduced for Manhat- 
ta n lattices by Gershensonl (|2005l ). and then generalized to arbitrary networks 



Ide Gier et al 



(|201lf ). In particular, the version of SOTL that we study is 
specifically designed to minimize the spatial heterogeneity of the density within 
the network. 

The remainder of this paper i s organized as follow s. In Section [21 we de- 



scribe the CA model introduced in Ide Gier et al.l ([20111 ). and define the network 



parameters we use for the simulations in this paper. Section [3] then defines the 
macroscopic quantities of interest in terms of observables of the CA model. In 
Section^ we describe the three traffic signal systems that we study in our sim- 
ulations, each of which has a different level of adaptivity. Then Sections [5] and [5] 
respectively describe the results of our simulations using time-independent and 



3 Sydney Coordinated Adaptive Traffic System 
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time-dependent boundary conditions. Section [7] discusses a modified version of 
the two-bin model in which the two bins are no longer assumed to be identical. 
Finally, Section [5] concludes with a discussion. 



2. Cellular Automata Model 

We briefly outline the cellular automata model used in our simulations, which 
we refer t o as the NetNaSch model. For a comprehensive description of the 
model see de Gier et al.l (|201 lh - 



Cellular automata (CA) are models which are discrete in time, space and 
state variables, whose dynamical rules are local. The NetNaSch model rep- 
resents a road network by a directed graph, in which the nodes represent in- 
tersections and the links represent streets. With each link is associated an 
ordered list of lanes, and each lane is a simple one-dimensional stochastic CA 
obe ying a (slight generalization of) t he Nagel-Schreckenberg (NaSch) dynamics; 



see 



Nagel and SchreckenbersJ (|1992T) . In addition, vehicles may move between 
neighboring lanes via simple lane-changing rules. Thus, the dynamics along 
each given link is essentially a standard CA freeway model, albeit with input 
and output rates that are determined dynamically by the rest of the network. 
The NetNaSch model intentionally avoids modeling the detailed motion of vehi- 
cles as they move through intersections; the underlying assumption being that 
the actual time a vehicle physically spends in an intersection is unimportant 
compared to the time spent on the inbound link waiting to traverse the inter- 
section. This course-grained approach allows the model to be easily applied to 
networks of arbitrary topology, using any choice of desired signal phasing. 

In order to mimic origin-destination behavior, the NetNaSch model demands 
that each vehicle makes a random decision about which link it wants to turn 
into at the approaching intersection. More precisely, for each node n, we assign 
to each ordered pair (1,1'), where / is an inlink and I' an outlink of n, the 
probability pr(l — > V) that a vehicle on / wants to turn into I' when it reaches 
n. The turning decision is made when the vehicle first enters I, since its choice 
of which link to turn into at the approaching intersection should influence its 
dynamics as it travels along In particular, it influences the vehicle's choice of 
when to change lanes. 

The NetNaSch model can be used with a variety of boundary conditions. In 
this paper we use open boundary conditions, and so the density in the network 
is not controlled directly. Instead, at each time step, vehicles enter and exit 
the network stochastically, according to prescribed input/output rates. We 
call in- and output at the boundary of the network exogenous, while internal 
sources/sinks are called endogenous, for example representing parking garages. 
In general, both exogenous and endogenous input and output are allowed. 

A boundary link is a link which has one of its two endpoints within the 
network, and one external to the network. Boundary links are classified as 
either boundary inlinks, if their to-node belongs to the network, or boundary 
outlinks, if their from-node belongs to the network. A bulk link is a link whose 
endpoints are both contained in the chosen network. In the NetNaSch model, 
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each lane A of each boundary inlink is assigned an input probability ay. at 
each discrete time step a new vehicle is inserted into the first cell of lane A 
with probability ax- Likewise, each lane A of each boundary outlink is assigned 
output probability j3\, which determines the probability that a vehicle wishing 
to exit the network from the last cell of lane A at a given time step actually be 
allowed to do so. 

The collections {a\}, {/3a} therefore specify the exogenous input/output of 
the network, i.e. they describe the level of demand imposed on the network by 
its environment. Intuitively, one can view a\ as being the density of an external 
reservoir of vehicles being fed into boundary lane A. And likewise, one can view 
(1 — /3a) as being the density of an external sink being fed vehicles by boundary 
out-lane A. Indeed, for the 1-lane NaSch freeway model, where the reservoir 
and sink can be thought of as on- and off-ramps, these interpretations are quite 
realistic. 

As discussed in Section [31 the boundary links are not considered to be part 
of our network, in the sense that we do not include their densities and flows in 
our network aggregated values. Instead, the boundary links are simply viewed 
as buffers allowing a realistic way to couple the bulk network to its external en- 
vironment. A practical issue that must be decided upon is how long to make the 
boundary links. There is no unique best answer to this question; a detailed dis- 
cussio n of the pros and cons of different possibilities is presented in lde Gier et al 



(|2011r ). For the simulations discussed in the present work, we simply set the 
length of the boundary links equal to the length of the bulk links. One advan- 
tage of this approach is that spillback caused by over-saturation on boundary 
outlinks is modelled dynamically, in a realistic way. 

In addition to these exogenous inputs/outputs, each lane of each bulk link is 
assigned an input probability 7a, and an output probability S\, which determine 
the rate of input and output from internal sources and sinks. For the simulations 
discussed in the present work, the internal sources and sinks were located near 
the middle of the lane, with the sink occurring before the source. 

For simplicity, we refer to the collection of all exogenous and endogenous 
inflow and outflow rates as the "boundary conditions" for the network, de- 
spite the fact that the endogenous rates are actually properties of the bulk. 
In principle then, the boundary conditions for a network are specified by the 
collections {a\ : A is a lane of an inlink}, {j3\ : A is a lane of an outlink}, {7A : 
A is a lane of a bulk link} and {^a : A is a lane of a bulk link}. For a given net- 
work, one could conceive of varying all of these parameters independently, from 
to 1, and studying the resulting distributions of flow and density. In order to 
meaningfully investigate MFDs however, we instead vary the a\, f3\, 7a and 8\ 
in a given systematic manner, corresponding to a reasonable demand scenario 
for an arterial network. We discuss several such scenarios in Section 12.41 We 
emphasize that the values of a\ 1 f3\, 7a and 5\ can vary with time. 

We now summarize the details of the specific network and input parameters 
simulated in the present study. 
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Figure 1: Illustration of a typical node in the simulated network. 



2.1. Links and lanes 

According to the NaSch model, the speed v of each vehicle can take one of 
« m „ + 1 allowed integer values v = 0,1,2,..., w max . Taking the length of a cell 
to be 7.5m, corresponding to the typical space occupied by each vehicle in a jam, 
and the duration of each time step to be Is, suggests w ma x = 3 is a reasonable 
choice for an urban network. I.e., each vehicle can move 0, 1, 2 or 3 cells per time 
step in such a CA model, depending on local traffic conditions. These are the 
values used in our simulations. In addition, the NaSch model (and consequently 
the NetNaSch model) includes, at each time step and for each vehicle, a random 
unit deceleration which is applied with probability p no ise- By setting p no i S c so 
that it is 0.5 when the current vehicle speed is w max , and 0.2 otherwise we obtain 
an average free-flow speed of approxi mately 60 km/hr, wh ich is typical of an 



arterial network. See the discussion in lde Gier et a l. (2011) for further details 



The particular network we simulated in this study consists of a regular 8x8 
square grid. Each link in the network has two lanes plus an additional right- 
turning lan^3. Fig. [TJ shows a typical intersection in detail. At the advice of the 
road authority in the Australian state of Victoria, the length of each bulk and 
boundary link was set to 750m, corresponding to 100 cells. This choice of link 
length corresponds to the distance between signalized intersections in an arterial 
network, and is typical of the suburban road network in Melbourne. The length 
of each turning lane was set to 120m. 



2.2. Phases 

Each node was given the same four phases: a north/south phase, an east/west 
turning phase, an east/west phase and a north/south turning phase. See Fig.[2j 
This fixed ordering of phases was applied to our simulations of SCATS. Note that 
the phase in Fig.[2]-(a) is not necessarily the first phase of the cycle; for SCATS 



4 Vchicles drive on the left side of the road in Australia. 
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(a) (b) (c) (d) 

Figure 2: The four phases used at each node of the simulated network. 



with signal linking this is determined by the linking protocol as described in 
Section PO 



2.3. Turning probabilities 

In all our simulations, each link was assigned the same turning probability of 
Pt for left and right turns, implying a probability l — 2p^ of continuing straight 
ahead. The probability p^ was set to 0.1 for the majority of our simulations. 
For comparison, in Section [5T3] we discuss simulations using = 0.15, 0.2. 

2.4- Boundary conditions 

We consider a number of representative scenarios for the boundary condi- 
tions, which we summarize as follows. 

I. Time-independent. All input and output rates are constant in time. Two 
main variations were studied. 

(a) Isotropic. The same value of a is applied to all inlinks, the same 
value of j3 to all outlinks, and the same values of 7 and S to all bulk 
links. To obtain MFDs, the values of a and j3 were varied while the 
values of 7 and d were held constant. 

(b) Anisotropic. The value of a on the west boundary is twice as large 
as on the other three boundaries, and likewise /3 is only half as large 
on the east boundary. This sets up a west-to-east anisotropy in the 
demand imposed on the network. The endogenous sources and sinks 
were ignored in this case, (7, S) = (0, 0). 

II. Time- dependent. In this case the boundary conditions are isotropic, but 
the values of a, (3, 7 and 8 are time-dependent. 

The values of a and /3 were changed every 30 minutes, and the network was 
simulated for 20 hours. For each of the three signal systems, the profiles 
of a and (3 were reverse-engineered to produce time series for the density 
pr ofiles that closely resemble the e mpirical data for Yokohama presented 
in iGeroliminis and Daganz 3 (120081) . In order to access the high density 



tails of the MFDs, howe ver, we allowed for slightly larg er values of the 
density than observed in Geroliminis and Daganzol ( 2008t) . Two extreme 



cases were studied. 
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(a) Boundary Loading. Vehicles can enter and exit the network only via 
boundary links. I.e. 7 and S are strictly zero. 

(b) Uniform Loading. We allow vehicles to be loaded and unloaded uni- 
formly throughout the network by setting "f = a and 5-/3. 

In the time-independent cases we compute MFDs by varying the exogenous 
demand, for a number of fixed choices of endogenous demand. This corre- 
sponds in some sense to viewing the endogenous demand as an inherent part 
of the network, comparable to the choice of signal system etc, and studying 
how the network responds to different levels of external demand. Systems with 
low endogenous demand correspond to "commuter corridors" ; portions of the 
arterial network wedged between a strong source (e.g. a residential area) and 
a strong sink (e.g. the central business district). Higher values of (7, 5) corre- 
spond to introducing shopping centers and parking garages etc into the bulk of 
the network. 

In the time-dependent case, in addition to the two extreme cases of boundary 
loading and uniform loading, we also studied a number of intermediate scenarios, 
in which the strength of the endogenous demand was non-zero but less than the 
exogenous demand. The resulting behavior of these systems was intermediate 
between the two extremes pi.a|) and (jll.bj) and so for the purposes of illustration 
it suffices to focus on bo undary and bulk loading. 

Finally, we note that iBuisson and Ladier (2009) present a careful discussion 



of possible sources of network heterogeneity in empirical studies of MFDs, in- 
cluding the position of detectors, types of roads, and traffic signals. We deliber- 
ately study symmetric square-lattice networks for which all links are equivalent 
and all nodes are equivalent. 

3. Observables 

We define the density, Pi(k), of link I at the kth time step of a simulation 
to be the fraction of all cells on I which are occupied at that instant. The flow, 
J\(k), of lane A during the kth time step is simply the indicator for the event 
that a vehicle crosses the boundary between a fixed pair of neighboring cells 
during the kth updatH The flow Ji(k) on link I at the fcth time step is then 
simply the sum of the J\(k) over all lanes A in link I. We emphasize that since 
our model is stochastic, the observables pi(k) and Ji{k) are random variables. 

Since we are interested in the dynamics on the order of traffic cycles, rather 
than iterations of our model, we bin the instantaneous link flow and density 
into bins of size b, using 6 = 5 minutes in our simulations. In a slight abuse of 
notation, we define 

kb kb 
fe=(t-l)6+l fe=(t-l)6+l 



3 In our simulations, the flow is measured 2v max cells from the upstream node of the link. 
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where the physical time t is measured in intervals of b = 5 minutes. 

Let us denote the set of all bulk links in the network by A. We emphasize 
that the boundary links in the NetNaSch model are not considered to be part 
of the network, and serve only as an effective means of connecting the bulk to 
the external environment. From the link observables (fTJ), we then define the 
following macroscopic network-aggregated observables: 



M*) ■= W|X|£M*)-p(*)] 2 , 

\ ^ (2) 

J W : =m£ J <W> 
M*) ^ £[•*(*) --W- 

Again, we emphasize that these macroscopic observables are random variables 
in our model, although by aggregating the data over both time and space the 
fluctuations of these macroscopic observables will be significantly suppressed 
relative to the original instantaneous link observables. 

The quantities J(t) and pit) are the network-aggregated flow and density. 
We refer to the quantities h p (t) and hj(t) as the heterogeneity (spatial vari- 
ability) of the density and flow respectively, since they give a measure of the 
extent to which the spatial distribution of the link-level observables differ from 
the corresponding network-aggregated values. Note that the heterogeneity lies 
strictly between and 1, and achieves the lower bound of only when the link 
observables are all equal. 

A fundamental question to be studied via our simulations is the extent to 
which pit) and/or h p (t) determine the value of J(t). The statement that an 
invariant MFD exists implies that Jjt) s h ould b e a f unction of pit) alone. How- 
ever, recent work by iMazloumian et aT ( 2010l ) and iGeroliminis and Sun TRB 



( 20111 ) suggest that h p (t) is also an important indicator of network performance. 
Our simulations confirm this. In fact, we find that both h P it) and hj(t) provide 
valuable indicators of network performance. 



3.1. Statistics 

For each distinct choice of traffic signal system and boundary conditions, 
we performed n independent simulations (with n ranging between 10 and 30), 
in order to estimate the expected values of the network- aggregated quantities 
defined in ([2]). For a given observable X(t), if we denote its realization in the 
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ith run by jW(t) then we compute 



err(X(i)) 



X(t) = -^X^(t), (3) 

(4) 



n '• , 
i=i 



a -7 rT -X(t)] 5 



where X(t) is the natural estimator for the expected value of X(t) and err(X(t)) 
is its standard error. 



4. Traffic Signal Systems 

We simulate and study the existence of MFDs in networks using three dis- 
tinct traffic signal systems: 

SCATS-L: A model of SCATS with linking and adaptive cycle lengths. 
SCATS-F: A "free" version of SCATS-L, with no signal linking. 
SOTL: Self-organizing traffic lights. 

The SCATS traffic signal system, which controls the traffic signals in nu- 
merous cities around the world, uses knowledge of the recent state of traffic 
to choose appropriate values of three key signal parameters: cycle length, split 
time, and linking offset. At each intersection it can adaptively adjust both the 
total cycle length, and the fraction (split) of the cycle given to each particular 
phase. In addition, it can coordinate (link) the traffic signals of several consec- 
utive nodes along a predetermined route by introducing fixed offsets between 
the starting times of specific phases, thereby creating a green wave. Both the 
SCATS-L and SCATS-F models are special cases of our general SCATS model, 
which we outline in Section 14.11 

As a benchmark with which to compare the SCATS-like traffic signal sys- 
tems, we als o considered the high ly-adaptive self- organizing traffic lights (SOTL) 
system (see de Gier et al.l 1 201 lh ) . The SOTL system is based on the simple 



principle that each node (intersection) should choose its current phase to be 
the phase which currently has the highest demand. Unlike SCATS, no direct 
coordination is enforced between the signals at neighboring nodes, however such 
coordination is often seen to arise via self-organization, since neighboring nodes 
do indirectly communicate with each other via the levels of traffic that they 
accept and release. The particular version of SOTL that we study here uses 
density as the demand metric, and it therefore strives to adaptively minimize 
the network's density heterogeneity. The details of the SOTL system are sum- 
marized in Section! 
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Figure 3: An 8 by 8 lattice network with 8 linearly linked subsystems. 

4.1. SCATS 

The SCATS control system adaptively controls three key signal parameters: 
linking offset, cycle length and split time. We discuss below how our model of 
SCATS chooses each of these parameters. 

Jj.,1.1. Linking 

A subsystem in a SCATS network is a group of nodes which all share a 
common cycle length. If a node does not belong to any subsystem, we call it 
a non-subsystem node. Within each subsystem, we appoint a unique master 
node m and a number of slave nodes s. Fig. [3] illustrates an 8 by 8 network 
with eight subsystems, each consisting of one master node and six slave nodes. 
To implement linking, each node is assigned a special phase V* , which is its 
linked phase. If the linked phase V* n of the master node m starts at time t then 
V* of the slave node starts at time t + T s , where T s is the offset. Ideally, the 
linking offset should be chosen based on the distance L between m and s, and 
the instantaneous local space-mean speed. In practice, actual implementations 
of SCATS tend to operate with fixed offsets during a given period of the day 
(for example morning peak hour). In our simulations we therefore use a fixed 
linking speed v = 54km/hr, which is just slightly less than the average free-flow 
speed of about 60km/hr, see Section [2TTI 

4- 1.2. Cycle length and split plan 

SCATS chooses the unique cycle length within a subsystem based on the 
local traffic conditions in the neighborhood of the master node, as quantified by 
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the Degree of Saturation (DoS). Every time a master node is about to restart its 
cycle, the cycle length is adjusted adaptively based on recent measured values 
of the DoS. In our model of SCATS, the cycle length is selected based on the 
volume ratio. For an inlink I and phase V the volume ratio is defined to be 

where V(l, V) is the measured traffic volume out of inlink I during phase T 3 , and 
S(V) is V's current split time. The quantity N(l, V) denotes a fixed benchmark 
volume for the link and phase, measured in vehicles per second^. If the volume 
ratio was large during the previous cycle, then the cycle length is increased by 
a fixed amount. Conversely, if green time was wasted during the previous cycle, 
the cycle length is decreased. The key underlying strategy is to attempt to 
keep the volume ratio within the range [0.85,0.95]. Once the cycle length is 
determined, the split of a phase is taken to be proportional to its traffic volume 
during the previous cycle. The specific details of cycle length and split time 
selection are discussed more fully in |Appcndix A.T 



4.1.3. Versions of SCATS 

Based on the model of SCATS outlined above, we considered two variants: 
SCATS-L and SCATS-F. The first variant, SCATS-L, operates on the linked 
network shown in Fig. [31 By contrast, SCATS-F operates on the network where 
no subsystems or linking are imposed, so that each node chooses its own cycle 
length and split time plan according to its local traffic state, independently of 
its neighbors. 

4.2. SOTL 

The third signal system we study is th e self-organizing traffic lights (SOTL) 



system described in Ide Gier et al.l (|201lh . While SCATS-L and SCATS-F are 



adaptive in their cycle length and split time selections, they both maintain a 
fixed cyclic ordering of each node's phases. By contrast, the SOTL system is 
acyclic, and is designed so that at each phase change, the phase which currently 
has highest demand is selected. This procedure is applied independently to each 
node. 

Suppose we agree on a suitable demand function d(V) which quantifies the 
demand of each phase V of each given node. Phases with large values of d(V) 
should be candidates for being the next choice of the active phase. However, one 
should also keep track of the time t(V) that each phase has been idle, since we 
do not want a given phase to remain idle for too long, unless it has strictly zero 
demand. The key idea behind SOTL is to compute a threshold function, k(V), 
for each phase P, which depends on both the phase's idle time and demand 
function, and when k{V) reaches a predetermined threshold value, 

k{V) > 6, 



In our simulations N(l,~P) was set to 1 veh/sec for all phases. 
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we consider making V t he active phase . For a detailed general discussion of the 
SOTL methodology, see lde Gier et all (|201lh . 

In the simulations performed in the current work, the demand d{V) of phase 
V was simply chosen to be the total number of vehicles over all its inlinks, and 
the threshold function was 



k{V) = 



d(V)r(P) 
E-P' d(V) 



(6) 



This particular choice for the demand function implies that SOTL attempts, at 
each instant of time, to adaptively minimize the network's density heterogeneity. 
We used a threshold value of 9 = 5. 

A precise algorithmic description of SOTL is given in Algorithm[2]in |Appcndix A. 2 



5. Simulations: Time-independent boundary conditions 

We simulated the network described in Section [5J using the three different 
traffic signal systems described in Section |H and measured the macroscopic 
observables defined in @. In this section, we present the results for the time- 
independent boundary conditions defined in Section 12.41 We simulated each 
system for 10 hours, which ensured that stationarity was always reached. 

5.1. Isotropic boundary conditions 

In this section, we present our results for simulations using the isotropic 
Boundary Condition (|I.a[) . In this case, we have two free parameters, a and 
/3, where a is the input probability on each boundary inlink, and (3 the output 
probability on each boundary outlink. Fixed values for the bulk input and out- 
put probabilities 7 and 5 were applied to all bulk links. Three different choices 
for the fixed values of (7, 6) were simulated: (7, 6) = (0, 0), (0.05, 0.1), (0.1, 0.2). 

We note that in the isotropic case, the SCATS-L system is somewhat ar- 
tificial, since there is no motivation for imposing linking if the demand is 
isotropic; we include SCATS-L here to enable comparisons with the results for 
the anisotropic network discussed in Section 15.21 

5.1.1. Zero 7 and 8 

Fixing 7 = 5 = 0, we simulated the network using a number of different 
values of a and /3, in order to obtain a range of values of the aggregated network 
density, ranging from very low to very high. We observed that, in all cases, the 
flow and density reach approximately stationary values by hours 5 or 6. For 
a given choice of traffic signals, the longest relaxation times were observed for 
flows close to capacity. For SOTL, stationarity was always achieved much faster 
than for SCATS (at comparable densities), typically by hours 3 or 4, which 
implies that the relaxation time of the network using SOTL is much smaller 
than when using SCATS. 

In Figs. [4(a)| [(b)] and [(cj] we plot J against p for SCATS-F, SCATS-L, and 



SOTL respectively, at hours 1, 2, . . . , 6 of the simulations. For each signal sys- 
tem, the low-density branch of the curves are essentially time-independent, as is 
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SCATS-F: isotropic and lime independent rates (y=0, 5=0), P T = 0.1 



SCATS-L: isotropic and time independent rates fy=Q, 5=0), p y = 0.1 
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Figure 4: Figs, [(a)] [(b)] and [(c)] show MFDs of SCATS-F, SCATS-L, and SOTL, at hours 
1,2, ...6, for a network with isotropic and time-independent boundary conditions, and no 
internal sources/sinks. Fig. |(d)| shows a comparison of the stationary MFDs for the three 
signal systems. Error bars corresponding to one standard deviation are shown, but are often 
smaller than the symbol size of the data point. 



the highly-oversaturated region of the congested branch. For intermediate val- 
ues of density, approximately in the range [0.3,0.8], we see a time dependence 
in the early hours of the simulation, however we also clearly see that the J vs ~p 
curves are converging to a well-defined stationary MFD as time increases. After 
approximately 6 hours of simulation, the J vs ~p curves at al l later times are 
essentially indistinguishable. We note that Mazloumian et al. ( 2010J ) observed 
similar time-dependent behavior during three-hour simulations of their model, 
and they concluded that such time-dependence would likely persist at all later 
times. Fig. [4] would suggest however that such time-dependence is in fact tran- 
sient, and, for all practical purposes, ceases to be o bservable after som e finite 
time. It is conceivable that the behavior observed in iMazloumian et al. ( 2010h 
is an consequence of their use of periodic boundary conditions. 
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SOTL: isotropic and lime independent rates (y=0, 5=0), p T = 0.1 

0.5 1 1 1 1 1 




0.2 0.4 0.6 0.8 1 

link density 

Figure 5: Single-link fundamental diagrams for a variety of links in the simulations of SOTL 
with isotropic and time-independent boundary conditions, at stationarity. The link labels are 
arbitrary, but a selection of representative links both near the boundary and well into the 
bulk have been chosen. Error bars corresponding to one standard deviation are shown. 



From Figs. |4(a)] and [4(b)] we can see that for both SCATS-F and SCATS- 
L, the instantaneous MFD curves appear to decrease monotonically with time 



towards their stationary limits. From Fig. 4(c) however, we notice that for 
SOTL the MFD at hour 1 lies below the limiting stationary curve. We shall 
return to a discussion of this transient behavior in Section 15 .1.31 

Fig. |4(d)| shows a comparison of the stationary MFDs for the three traffic 
signal systems, SCATS-F, SCATS-L, and SOTL. In the low density regime, the 
performance of each system is quite similar. However, SOTL clearly allows the 
network to reach a higher capacity: SOTL achieves a maximum flow 0.428 ± 
0.001 at a density around 0.341 ± 0.001, while SCATS-F and SCATS-L obtain 
maxima of 0.390 ±0.001 and 0.388 ±0.004 at densities 0.184 ±0.004 and 0.187 ± 
0.009, respectively. This represents a 9.6% increase in network capacity by using 
SOTL, compared with SCATS. We note that SCATS-F performs better than 
SCATS-L in the high density regime. This is in fact to be expected, since for 
an isotropic network, linking should at best be merely unhelpful, and at worst 
it will be counterproductive because it will reduce the system's adaptivity. We 
return to the comparison between SCATS-F and SCATS-L in Section [OJ 

To facilitate comparison with empirical data, we note that, using the link 
and cell lengths described in Section [2~T1 to convert the above density and flow 
estimates to more standard units, one multiplies the quoted densities by 75 to 
obtain densities in units of vehicles per kilometer, and multiplies flows by 1800 
to obtain flows in vehicles per hour. E.g., the SCATS-F values above correspond 
to a maximum flow of 698 ± 7 veh/hr at a density of 14.3 ± 0.7 veh/km , which 
is c omparable to the empirica l valu es reported in iBuisson and Ladier 
and iGeroliminis and Daganzo (l2008lh Analogous conversions can be applied 
throughout. 
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(d) Flow heterogeneity: SOTL 



Figure 6: Heterogeneity versus density at hours 1,2,. . . ,6 for network governed by SCATS-L 
and SOTL, with isotropic and time-independent boundary conditions. Error bars correspond- 
ing to one standard deviation are shown, but are often smaller than the symbol size of the 
data point. 



5.1.2. Comparing individual link fundamental diagrams 

For comparison with the MFDs shown in Fig. |U in Fig. [5] we plot the 
FDs for a number of representative links in a network governed by SOTL, 
at stationarity. It is known that the NaSch model with time-independent 
bounda ry conditions gives a triang ular fundamental diagram at stationarity; 
see e.g. ISchadschneider et al. ( 2011 ). However, for our network model the links 
are essentially correlated NaSch models with time-dependent and random a and 
f3, so it is not clea r whether the lin ks should have FDs, and if so, what form 
they should have. IWu et al.l (|2011l ) recently performed an empirical study of 
"arterial fundamental diagrams" , i.e. fundamental diagrams for single links in 
arterial networks, and argued that the presence of traffic signals should imply 
that the usual triangular FD be replaced with a trapezoidal FD. The results 
presented in Fig. [5] confirm this. 
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We note that there is clearly significant scatter in the link FDs in Fig. [5j 
emphasizing that even networks with homogeneous infrastructure and isotropic 
demand can have inhomogeneous spatial distributions of traffic. 

5.1.3. Heterogeneity 

In order to gain insight into the underlying cause of the transient behav- 
ior displayed by the MFDs, we consider how the corresponding heterogeneity 
curves evolve with time. We begin by considering th e density heterogeneity. In 



Fig. [6(a) we plot h p for SCATS-L at hours 1,2,. . .6. iGeroliminis and Sun TRB 



(2011) find empirically that for a given network density, the flow should be 



lower when the density heterogeneity is higher. Comparison of Fig. |4(b")1 with 



Fig. 6(a) confirms this; if one observes the time evolutions of the MFD curve 
and hp vs p curve in the neighborhood of a given value of p, there is a clear 
anticorrelation between the flow and the density heterogeneity. At moderate 
values of the density, h p starts at low values early in the simulation, and then 
increases to a stationary curve at around hour 5 or 6. Conversely, the flow starts 
at relatively high values and decreases to its stationary value, again reaching 
stationarity at around hour 5 or 6. Similar behavior is observed for SCATS-F. 



If one considers the analogous plots for SOTL, Figs. 4(c) and 6(c) at hours 
2,. . . , 6, precisely the same behavior can be observed, although it is less pro- 
nounced. The most obvious feature in Figs. 4(c) and 6(c)[ however, is the 



behavior of the instantaneous curves at hour 1: for moderate values of density, 
the flow is below its stationary limit, while h p is above its stationary limit. 
This again demonstrates the anticorrelation between the network flow and the 
density heterogeneity. In fact, the same type of behavior displayed by SOTL at 
hour 1 is also present in SCATS-L and SCATS-F, however it occurs before hour 
1 in these cases and so is not visible in Figs. 4(a)| |4(b)| and 6(a) This type of 



behavior is presumably related to the fact that the simulations are started from 
a completely empty network, with vehicles only entering via the boundary. 

The general conclusion to be drawn from these observations is that while 
the transient behavior of the flow at a given density might be subtle (e.g. non- 
monotonic in time), at all times there is a strong anticorrelation between flow 
and density heterogeneity. We shall see in Section \5\ that this relationship be- 
tween the transient behavior of J and h p plays an important role in understand- 
ing hysteresis in networks with time-dependent boundary conditions. 

As a practical observation, we note that h p is considerably smaller when 
the network is governed by SOTL than when governed by SCATS-L, which is 
unsurprising given that SOTL is specifically designed to adaptively reduce the 
density heterogeneity. Given the previous observation that the SOTL MFD lies 
strictly above the SCATS-L MFD, this suggests that a methodology of adaptive 
density homogenization can provide a highly effective means of network control. 

The relationship between the transient behaviors of the flow and the flow 
heterogeneity is more complicated, however. We shall present a careful time 
series analysis of the cross correlations between hj and J elsewhere. However, 
the stationary behavior of the hj vs ~p curves appear to contain a considerable 
amount of information. In particular, for all three signal systems studied, there 
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appears to be a point of inflection in the stationary hj vs p curve either at, or 
slightly before, the density p c at which the network reaches capacity. Figs. |6(b)| 



and 6(d) show hj vs p for SCATS-L and SOTL, respectively. For the case of 
SOTL there is in fact a long plateau, implying that for the majority of the free- 
flow regime the flow heterogeneity is independent of the network density in this 
case. 

5.1.4- Effects of internal sources/sinks 

We now consider the effect of introducing non-zero values of the internal 
input and output probabilities, 7 and 6, which model endogenous sources/sinks 
such as parking garages. Figs. [7(a)1 and [7(b)] show the MFDs for SCATS-L and 



SOTL produced by varying a and /3 for two distinct, non-zero, fixed values of 
(7,(5), and compares them with the (7, 5) — (0,0) case already discussed. The 
results for SCATS-F arc intermediate between those of SCATS-L and SOTL. 
The first observation is simply that different values of (7, S) clearly produce 
different MFD curves as a and /3 are varied. This is to be expected, since the 
stronger the internal sources and sinks, the more homogeneous is the network; 
for sufficiently strong internal sources and sinks the effects of the boundary 
become essentially irrelevant. 

The second observation is that for both SCATS-L and SOTL, the MFDs are 
translated to the right (higher densities) as 7 and 5 increase. This is intuitively 
reasonable; for (7, S) strictly zero the links deep in the bulk of the network would 
be expected to have lower density than those near the boundary, implying that 
the aggregated network density should be lower. This scenario explains why 
the translation of the MFD is more pronounced for SCATS-L than for SOTL, 
since we have already seen in Fig. |6] that the heterogeneity of SCATS-L is much 
higher than that for SOTL. 

A final observation is that for SOTL the capacity is also marginally higher 
with non-zero 7, 5 than with zero 7, 5. For SCATS-L, there is considerably more 
statistical noise in the high density branch when simulating with non-zero j,5, 
and it is not clear to what extent the capacities change, if at all. 

Apart from these translations and rescalings, Figs. 7(a) and |7(b)j show that 



no qualitative change in the shapes of the stationary MFDs is introduced by 
applying non-zero internal sources/sinks. By contrast, we will see in Section [6] 
that when using time-dependent boundary conditions, the introduction of suf- 
ficiently strong internal sources and sinks can qualitatively affect the behavior 
of the network. 

5.1.5. Turning probabilities 

We now briefly turn our attention to the impact on MFDs of varying the 
turning probabilities, again using isotropic boundary conditions (Boundary Con- 
dition dLi|). Figs. |8(a)| |8(b)| and |(cj compare the MFDs produced for net- 



works with turning probabilities p T = 0.1,0.15,0.2, for SCATS-F, SCATS-L, 
and SOTL, respectively. While the low density branches are invariant with p^, 
it is clear that the high density branches become more rapidly decaying as px 
increases. For both pj- = 0.15 and pr = 0.2, there is a critical density pj such 
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Figure 7: MFDs for SCATS-L and SOTL with 7 = 0, 0.05, 0.1 and 5 = 0, 0.1, 0.2. Error bars 
corresponding to one standard deviation are shown. 



that the flow remains constant for all ~p > p y This value of pj decreases as px 
increases. Although it is not observable for px = 0.1, it is presumably present in 
principle but at a density higher than any we simulated. We note that even for 
p > /Oj, the flow is non-zero, so the network is not locked into a rigid gridlock. 
The value of this plateaued high-density flow can be seen to decrease as px 
increases. This is all in accordance with one's intuition that the network should 
behave less efficiently as the probability for drivers to make turns increases. 

On a practical level, we note that for px = 0.15 and px = 0.2 the capacities 
for SCATS-L and SCATS-F are again very close, and again both are significantly 
lower than the corresponding capacity observed for SOTL. 

Fig. |8(d)| shows the instantaneous MFDs at hours 1,2,. . . ,6 for a network 
governed by SCATS-F, with px = 0.2. The transient behavior is qualitatively 



the same as that shown in Fig. 4(a) for px = 0.1. The low density branch 
is stationary already by hour 1, while the high density branch decreases to its 
stationary limit, which is reached by hour 4 or 5. The extent to which the curves 
at hours 1 and 2 differ from their stationary limits is clearly much larger for 
px = 0.2 than observed for px = 0.1 however. Similar behavior was observed 
for SCATS-L and SOTL. 



5.2. Anisotropic boundary conditions 

In this section we present our results for simulations using Boundary Con- 
dition (jl.bjl . In this case, we again have two free parameters, a and /3, where a 
is the input probability on each inlink on the west boundary, and (3 the output 
probability on each outlink on the east boundary. All other boundary inlinks 
have input probability a/2 and all other boundary outlinks have output proba- 
bility 2/3. No internal sources or sinks are present (7 = 6 — 0) and the turning 
probability is px = 0.1. These boundary conditions imply that the demand in 
the west-to-east direction is twice that of other directions. In the presence of 
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SCATS-F: isotropic and time independent rates (y=0, 5=0), time 6hr 



SCATS-L: isotropic and time independent rates fy=0, 5=0), time 6hr 
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Figure 8: Figs. [(a)] [(b)] and[(c)]show MFDs of SCATS-F, SCATS-L, and SOTL, with turning 
probability px = 0.1,0.15,0.2, for a network with isotropic and time-independent boundary 
conditions. Fig. |(d)| shows MFDs for SCATS-F at hours 1,2,. . . ,6 with p T = 0.2. Error bars 
corresponding to one standard deviation are shown, but are often smaller than the symbol 
size of the data point. 



such anisotropy, applying linking with SCATS-L is a very natural thing to do, 
and our simulations of SCATS-L are linked in the west-to-east direction. 



In Fig. 9(a) we plot a comparison of the stationary MFDs of the networks 
using SCATS-F, SCATS-L, and SOTL, corresponding to hour 6 of our simula- 
tions, by which time each system had reached approximate stationarity. The 
transient behavior prior to stationarity is qualitatively the same as described for 
the case of isotropic boundary conditions discussed in Section I5TT1 The first ob- 
servation is simply that for each signal system, well-defined stationary MFDs do 
exist in the presence of anisotropic boundary conditions. However, the shapes 
of the stationary curves are clearly quite different to those produced by isotropic 
boundary conditions as shown in Fig. 01 While in the isotropic case the MFDs 
are smooth curves which qualitatively resemble a typical single-link FD, the 
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anisotropic and time independent rates (y=0, 5=0) , p T = 0.1 at 6 hr 
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Figure 9: Comparison of SCATS-L, SCATS-F, and SOTLat stationarity, on a network with 
anisotropic a and in the west-to-east direction. Linking in SCATS-L is applied along the 
high-demand direction. Fig. |(b) [ shows the MFDs while Fig. |9(b)| shows density heterogeneity. 
Error bars corresponding to one standard deviation are shown, but are often smaller than the 
symbol size of the data point. 



MFDs shown in Fig. 9(a) have a more intricate structure. In particular, for 
all three signal systems, the stationary MFD displays two steep drops in flow. 
The first of these occurs for p w 0.4, and is most pronounced for SCATS. For 
SOTL, the location of this first drop is close to the density p c which produces 
maximum flow, while for SCATS it occurs well beyond p c . 

As observed in the isotropic case for pt = 0.2, for densities above a critical 
value, pj, an oversaturated regime with constant flow is obtained. Again, even 
for ~p > pj the flow remains non-zero, although it is very low. For all three 
signal systems, a second steep drop can be observed for ~p « p y This drop is 
very sharp, even for SOTL. Some insight into this behavior can be obtained 
by comparing with the corresponding density heterogeneities in Fig. |9(b)| For 
all three signal systems, the location of the drop in the MFD near pj coincides 
exactly with a cusp in the density heterogeneity. In addition, for SCATS-L and 
SCATS-F, the sharp drop in the MFD near ~p » 0.4 coincides exactly with a 
sharp rise in the density heterogeneity. In summary, we see that anisotropies in 
demand significantly affect the shape, but not the existence of the MFD curves. 

We now turn to a comparison of the different signal systems. We begin 
by comparing SCATS-L and SCATS-F. For networks with strongly anisotropic 
demand, one would intuitively expect signal linking to play a beneficial role. 
Indeed, if one focuses only on travel times along the peak (linked) direction, 
the SCATS-L system performs considerably better than SCATS-F in the low 
density regime. However, we see from Fig. 9(a) that the MFDs produced by 
SCATS-L and SCATS-F are essentially indistinguishable. Although SCATS-L 
provides improved performance for the links in the peak direction, this benefi- 
cial treatment of the linked subsystem occurs at the expense of the non-linked 
routes, and no net benefit is observable when using the network MFD as the 
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performance metric. At high densities, we find no significant difference between 
either the network MFDs or the travel times along the linked direction when 
using SCATS-L or SCATS-F, in contrast to the isotropic case (c.f. Fig. |4(d)[ ), 
where linking at high densities was clearly counterproductive with respect to 
the MFD. 

Finally, let us compare the performance of SOTL with that of the SCATS 
systems. At densities lower than 0.2 and higher than 0.4 all three systems per- 
form similarly, although we note that at low densities SOTL has both higher 
network flows and lower travel times in the linked direction than either SCATS- 
L or SCATS-F. For moderate densities, however, we observe that SOTL has 
significantly higher network flow, and lower density heterogeneity, than cither 
SCATS system. This suggests that a methodology of adaptive density homog- 
enization can provide a more effective means of network control than signal 
linking, evne in in networks with inherently anisotropic demand. 



6. Simulations: Time-dependent boundary conditions 

In the previous section, we applied time-independent boundary conditions, 
and observed that the system took up to 6 hours to reach stationarity. In 
real-world scenarios, however, traffic demand typically varies with time, and 
in practice a network may never actually reach stationarity. Understanding 
transient behavior of the network dynamics is therefore of significant practical 
importance. 

In this section we present our results for simulations using Boundary Con- 
dition [IT] of Section 12.41 We select appropriate values of a and (3 so that we 
can simulate the network traffic variation during a typical weekday. Specifically, 
we consider a 20 hour period, and enforce two peaks in the demand; one cor- 
responding to the morning, the other to the afternoon. In order to access the 
high density regime of the fundamental diagram, the average network density 
in the afternoon peak is selected to be around 0.6 to 0.7. At each instant, the 
same value of a (resp. /3) is applied to each boundary inlink (resp. outlink), so 
the boundary conditions are isotropic. We consider two scenarios for the inter- 
nal sources and sinks: 7 = 5 = 0, in which case the demand is entirely driven 
by the boundary; and 7 = a, 8 = fi, in which case the bulk input (output) 
occurs at the same rate as the boundary input (output). We refer to these two 
cases as boundary loading and uniform loading, respectively. Unlike the time- 
independent case discussed in the previous section, in the time-dependent case 
we find that varying the relative strength of the internal sources/sinks compared 
with the boundary demand can produce rather different qualitative behavior. 

The time series of the resulting network density and flow for SCATS-L are 
shown in Fig. [TO] the corresponding profiles for SCATS-F and SOTL were very 
similar, although SOTL produced higher and sharper peaks. For each of the 
three signal systems, and for both boundary and uniform loading, the density 
and flo w profiles were engineered to clo sely resemble those of the Yokohama 



study in Geroliminis and Daganzo (|2008t ) , however we allowed for slightly larger 
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Figure 10: Time series of network-aggregated flow and density for the SCATS-L traffic signal 
system under time-dependent isotropic boundary conditions. Error bars corresponding to one 
standard deviation are shown. 



values of the density in order to be able to access the high density tail of the 

MFD. 

Figs. ll(a)| H( c ) an d H( e ) show the relationship between the density and 
flow for each of the three signal control systems studied, in the case of uniform 
loading, while Figs. 12(a)| 12(c) and 12(c)| show the analogous plots in the case 
of boundary loading. Each plot is in fact a parametric plot in the density-flow 
plane, parameterized by time. In each case, the density-flow curve obtained dur- 
ing the first 4 hours of the simulation coincides exactly with the corresponding 
stationary MFD. During this period, the network is initially empty, and as the 
density increases the flow increases until maximum flow is obtained. Through- 
out this period, the network remains uncongested, and no transient effects are 
observed. However, as the morning peak in demand is approached, which occurs 
at around hour 6 of the simulations, we see that the density-flow curves develop 
non-trivial time dependences, and hysteresis effects emerge. 



6. 1 . Hysteresis 

Let us analyze the observed hysteresis patterns in more detail. We begin with 



the case of uniform loading. Fig. 11(e) shows the time evolution of the MFD for 
SOTL. The system reaches capacity at approximately hour 4, after which time 
flow decreases as density is further increased. This continues until around hour 
6, when the first peak in demand is reached. After this point, density drops 
and flow begins to increase again until it again reaches capacity, at around 
hour 8. However the curve followed during this "recovery" process (hours 6 
to 8) lies below the curve followed in the original "loading" process (hours 4 
to 6). The MFD curve from hour 4 to hour 8 therefore defines a clockwise 
hysteresis loop. Similar behavior occurs at later times also. Indeed, we observe 
a further two clockwise hysteresis loops: from capacity, to low density, back to 
capacity (hours 8 to 12); from capacity, to high density, back to capacity (hours 
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Figure 11: Performance of network under time-dependent isotropic boundary conditions, when 
using SCATS-F, SCATS-L, and SOTL traffic signal systems with uniform loading. Left col- 
umn: Instantaneous MFDs. Right column: Density heterogeneities. Error bars corresponding 
to one standard deviation are shown. 
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SCATS-F: time dependent rates, p t = 0.1 (7 = 0, 5 = 0) 
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Figure 12: Performance of network under time-dependent isotropic boundary conditions, when 
using SCATS-F, SCATS-L, and SOTL traffic signal systems with boundary loading. Left col- 
umn: Instantaneous MFDs. Right column: Density heterogeneities. Error bars corresponding 
to one standard deviation are shown. 
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12 to 16). Finally, the system recovers from capacity to the very low density 
regime (hours 16 to 20); this final recovery curve is initially below the initial 
loading curve (hours to 4), but coincides with it at sufficiently low densities. 
Qualitatively similar behavior is also observed for SCATS-L and SCATS-F, 
although in these cases the system during recovery generally does not actually 
return to the original loading curve, so that although hysteresis is clearly present, 
closed loops are not necessarily formed; this is particularly true for SCATS-F. 

Clockwise hysteresis loops we re recently observed in an empirical study of 
MFDs in Toulouse, presented in iBuisson and Ladier (l2009h . where it was ar- 
gued that hysteresis is caused by spatial heterogeneity in network density. In 
Figs. 11(b) 11(d) and 11(f) we show plots of the density heterogeneity, under 
uniform loading, for each of the three signal systems studied. In each case, the 
heterogeneity displays hysteretic behavior very similar to that observed in the 
corresponding MFD, but with opposite orientation; the MFDs in Fig. [11] all 
display clockwise hysteresis loops, while their corresponding heterogeneity plots 
display anticlockwise loops. This behavior is to be expected when loading is 
uniform; the level of density heterogeneity produced when a network is loading 
is governed predominantly by the locations of the sources, whereas the level of 
heterogeneity produced when the network recovers will depend not only on the 
locations of the sinks, but also on drivers' behavior as they disperse through 
the network towards those sinks. Therefore, whenever the spatial distribution 
of the sources is sufficiently uniform, one would expect that the heterogene- 
ity of loading should be s maller than that of re c overy . A recent study of the 
two-bin model presented in iGavah and Daganzol (|201ll ). under the assumption 
of perfectly symmetric loading, provides a simple analytical explanation of this 
behavior. 

If the constraint of perfectly uniform loading is relaxed however, different 
behaviour can arise. If loading is sufficiently non-uniform then density may in 
fact become more heterogeneously distributed in loading than in recovery, which 
would produce anticlockwise hysteresis loops in the MFD, rather than clockwise 
loops. To study this possibility, we have therefore simulated networks in which 
the boundary sources/sinks are stronger than the internal sources/sinks. Situ- 
ations where such behavior might arise include commuter corridors as well as 
arterial networks in the presence of perimeter control, or gating. The most ex- 
treme example of such loading is when the strength of the internal sources and 
sinks is set identically to zero, corresponding to the profiles in FigfT2l 

Consider Fig. 12(c) for the SOTL signal system, where the effects are most 
pronounced. During the first 10 hours of the simulation corresponding to the 
lead-up to, and recovery from, the morning peak hour, we observe two distinct 
hysteresis loops. The first of these loops, which occurs at moderately high den- 
sity, is traced out by the flow-density curve in an anticlockwise direction, as 
time evolves, while the second loop, occuring at low density, is traced out in a 
clockwise direction. As the density then increases again due to the afternoon 
peak, a second, larger, anticlockwise hysteresis loop is traced out at high densi- 
ties. The behavior at moderately high density can be readily understood from 
Fig. 12(f) as follows. An initially empty network is loaded just above its ca- 
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pacity, in a non-uniform manner which leaves some links very congested, with 
others still uncongested. After the morning peak (around hour 6), as the in- 
flow drops and outflow increases, vehicles disperse through the network, and the 
density becomes more evenly distributed. This causes the anticlockwise loop at 
moderately high density as well as the larger anticlockwise loop at high den- 
sity. In Section [7] we provide a simple explanation for the se observations using 
a hete rogeneous version of the two-bin model discussed in Gavah and Daganzol 



(2011). 



The behavior of Fig. [T2(f)| at low density requires some further consideration. 



From Fig. 12(f) we can see that instead of the density heterogeneity displaying 
an anticlockwise loop during the low density clockwise hysteresis loop of the 
corresponding MFD (hours 8 to 12), it displays a figure-eight pattern. (A similar 
effect can be seen shortly after hour 16; an initially sharp drop in heterogeneity 
followed by an increase.) This at first seems counterintuitive; around p « 0.3, 
both J and h p are lower during recovery than during loading. However, h p 
only quantifies the extent to which densities vary between links; it tells us 
nothing about how the density varies along the individual links. As the density 
decreases from capacity, the downstream queues on the congested links near the 
boundary start to break up as traffic flows out of the network; this can have 
a significant effect on the density heterogeneity, without necessarily affecting 
the transient values of the flows, which are measured at the upstream end of 
the links where the congestion is not yet changing considerably. Eventually 
however, this transient effect abates and the recovery heterogeneity curve rises 
above the loading heterogeneity curve. 

The two scenarios, boundary loading and uniform loading, can be considered 
opposite extremes; in practice, one would likely expect that the strength of the 
internal sources and sinks would be non-zero but smaller than the boundary 
rates. We therefore also simulated a scenario with 7 = a/3 and 6 = /3/3. The 
results are intermediate between the uniform loading and boundary loading 
cases discussed above. In particular, for the SOTL system, we observe both 
anticlockwise and clockwise hysteresis at high density: anticlockwise hysteresis 
during hours 4 to 8 (as observed for boundary loading) ; and clockwise hysteresis 
during hours 12 to 16 (as observed for uniform loading). 

Finally, let us compare the behavior of the SOTL simulations with the corre- 
sponding simulations of SCATS. For both boundary loading and uniform load- 
ing, we again find that SOTL has lower values of density heterogeneity and 
higher capacities than both the SCATS systems. For the specific case of bound- 



ary loading, we note that while Figs. 12(a) and 12(c) display clockwise loops in 



the MFDs of SCATS-L and SCATS-F, as observed for SOTL, the anticlockwise 
loops appear to be absent. There is, in fact, a small, but well-defined, clockwise 
hysteresis loop in the density heterogeneity plot for SCATS-L, and if one zooms 



in on Fig. 12(c) then a corresponding tiny anticlockwise loop can be observed 
in the MFD. However the loop is smaller than the size of the error bars in the 
simulations and so its existence cannot be firmly established. The presence of 
strong anticlockwise loops for SOTL but not for SCATS can be understood as 
a consequence of SOTL's generally stronger ability to homogenize the network 
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density. 



7. Two-Bin model 



The two-bin model was introduced by iDaganzo et. ~a\ ( 2011 ). It represents 



the interaction between two subnetworks, or bins, in a larger road network. 
A state of the system consists of the pair {pi,p-z) : and the expression for the 
flow J(-) in each bin is assumed to be given by the same triangular MFD, with 
capacity (p c , J c ) and jam density p y The dynamical evolution of the system is 
defined by the following system of ordinary differential equations 



d/9l _ Ql - biJ(pi) +P2J{P2) -piJ(pi) 

dt Lx 
dpi _ a.2 - b 2 J(p 2 ) +p\J{pi) ~P2J{P2) 
dt L 2 



(7a) 



for pi , p 2 < Pj , and 

dpi dp 2 

-dr = ^r = 0lipiOiP2 = p >- (7b) 

The model as stated in (J7]) has 8 free parameters: Oj, the rate of inflow into 
bin i; bi, the proportion of traffic flowing out of the network from bin i; Pi, the 
proportion of traffic turning out of bin i into the other bin; Lj, the total network 
length of bin i. 

The perfectly symmetric case, in which a\ = 02, by = 62 , pi = P2, and 
Li = L2, was studied in detail by iGavah and Daganzo ( 2011 ). It was found 



that if hysteresis was observed in the aggregated MFD, it was typically oriented 
clockwise. The results for the uniform loading scenario studied in Section [5] are 
in perfect agreement with these results. 

In this section, we consider instead a highly asymmetric case, in which a 2 = 
62 = 0, and pi and P2 are not necessarily equal. This can be viewed as a 
two-bin model of the boundary loading scenario studied in Section [6] In this 
interpretation, bin 1 corresponds to the links adjacent to the boundary, while 
bin 2 corresponds to the remainder of the bulk links, see Fig. [3] This implies 
that L\ is slightly smaller than L2. The aggregated density and flow are then 

p _ Lipi + L2P2 j _ L l J{p 1 ) + L 2 J(p 2 ) ^ 
L\ + L2 Li + L 2 



Figs. 13(a) and |13(b)] show phase plots for the system (J7| during loading and 



recovery, respectively, for a typical choice of the model parameters. The loading 
process has ai — 0.1 and 61 = 0, while the recovery process has ai = and 
bi = 0.1, so that no vehicles exit the network during loading, and no vehicles 
enter the network during recovery. In b oth cases pi = 0.08, P2 = 0-02. For 
the perfectly symmetric case studied by IGavah and Daganzol ( 2011 ). the line 



Pi = P2 corresponds to a line of unstable fixed points, however we note that for 
the asymmetric system studied here perfectly balanced states (with p\ = p 2 ) 
never occur at late times. 



2N 



(a) Loading: 01 = 0.1, f>i = 0, pi = 0.08. (b) Recovery: ai = 0, h = 0.1, pi = 0.08. 



a, = 0.1, b 1 =0.1, p, = 0.08, p ? = 0.02 a, =0.1, b, = 0, p, =0.12, p 2 = 0.02 




(c) MFD corresponding to |(a)| and |(b)| (d) Loading: ai = 0.1, fei = 0, pi = 0.12. 



Figure 13: Solutions to (0 with ci2 = = 0. The loading paths have fri = while the recovery 
paths have a\ = 0. Fig. |(c)| shows the trajectory in the p — J plane corresponding to a typical 
loading/recovery process, given by the black paths in |(a)| and |(b)| The star represents the 
maximum value of density obtained during this loading/recovery process, and the diamond 
(triangle) represents the maximum flow attained during loading (recovery). Fig. |(d)| shows 
the effect on Fig. |(a)| of increasing p\. 



To illustrate the effect on the MFD of combining such loading and recov- 
ery processes, let us now consider the black trajectories shown in Figs 
and p[b)j In Fig. \HjsT 



13(a) 



0, 



the black trajectory starts with p\ > and pi 
which corresponds to the state of the boundary loading scenario simulated in 
Section [S] during the early stages of the simulation. Suppose we stop the load- 
ing process at a finite time, denoted by the star in Fig. 13(a)| and then be- 
gin a recovery process. The resulting recovery trajectory is shown in black in 
Fig. 13(b) The diamond (triangle) in Fig. 13(a) (Fig. 13(b) ) shows the location 
at which maximum flow was attained on the black loading (recovery) trajectory. 
In Fig. 13(c) we show the trajectory in the p— J plane resulting from combining 
these loading and recovery processes, superimposed on the underlying MFD of 
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the model. There is a clearly visible anticlockwise hysteresis loop as the sys- 
tem recovers from high density, precisely as observed for the boundary loading 
scenario simulated in Section [5] 

Since the recovery path in Fig. |13(b)| initially moves towards more balanced 
states, i.e. moves toward the diagonal, the higher flow during the initial stages 
of recovery can be seen as a consequence of lower heterogeneity, as observed 
in Section [6] Note that the location of the maximum flow during recovery 
occurs precisely at zero heterogeneity, p\ — pi. As the recovery trajectory 
continues further, Fig. |13(b)| shows that the heterogeneity increases again, and 
the corresponding trajectory in the MFD then drops below the original loading 
curve and a clockwise hysteresis loop results. This combination of anticlockwise 
hysteresis loops at high density and clockwise loops at low density agrees exactly 



with the behavior observed in Fig. 12(e) This behavior is quite generic, and 



similar MFD hysteresis is observed for any similar pairs of loading and recovery 
trajectories, whenever loading begins with pi ^> p 2 and ends with pi ^> p c . 

As a final observation, Fig. |13(d)| shows an alternative loading process, for 
which all parameters are the same as Fig. 13(a) except that p\ = 0.12 is slightly 
higher. While the trajectories are qualitatively similar, it is apparent that even 
a small variation in the rate that vehicles from the boundary can enter the 
network can produce quite different loading trajectories, starting from a given 
initial state. One significant factor that would contribute to the effective value 
of pi in an actual arterial network is the type of traffic signal system used. From 
the perspective of the two-bin model, it is therefore unsurprising that different 
types of signal systems can have rather different levels of hysteresis in their 
MFDs. 



8. Discussion 

We have studied macroscopic fundamental diagrams (MFD) on a square- 
lattice traffic network for a variety of traffic scenarios. In particular we have 
studied networks as a function of both anisotropy (directional bias) and unifor- 
mity (presence of bulk sinks and sources) in demand. We furthermore studied 
the impact of various turning probabilities, and we have studied these networks 
with three different traffic signal systems, SCATS-F, SCATS-L and SOTL. 

Our main findings include the following: 

• For time-independent demand, MFDs do exist even when demand is not 
uniform, but their shapes depend on the nature of the non-uniformity: 

— Systems with more uniformly distributed demand achieve similar ca- 
pacities, but capacities occur at higher densities; 

— Systems subject to anisotropic exogenous demand display a steep 
drop in the flow just beyond the maximum of the MFD; 

— As turning rates increase, capacity and jamming occur at lower den- 
sities, capacity decreases, and the steepness of decay in the congested 
branch of the MFD increases; 
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• For time-dependent demand, MFDs show clear hysteresis which is strongly 
correlated with the spatial heterogeneity of the density. The qualitative 
behaviour of this hysteresis is strongly dependent on the level of uniformity 
of loading and unloading; 

• The choice of traffic signal system plays a crucial role in determining a 
network's performance. The idealized control scheme SOTL, which is 
designed to uniformize the network density distribution, always results in 
a higher MFD compared to the commonly used SCATS system. SOTL 
increases the network capacity and produces higher flows in the congested 
regime. 

We finally remark that given the strong dependence of the choice of traffic 
signal system on the shape of the MFD, one can in principle use MFDs as a 
metric for comparing the performance of different traffic signal systems. 
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Case 1: 


if 


Case 2: 


if 


Case 3: 


if 


Case 4: 


if 


Otherwise: 



Appendix A. Details of Traffic Signal Control Systems 

Appendix A.l. SCATS 

The strategy for adapting the cycle length C based on the volume ratio R, 
defined in ([5]), of a master node is presented as below. 

Algorithm 1. SCATS cycle length decision 

C = MIN &R> 0.4, then C = STOPPER 
C = STOPPER &R< 0.2, then C = MIN 
R > 0.95, then C = min{C+ STEP, MAX} 
R < 0.85, then C = max{C- STEP, STOPPER} 
C remains unchanged. 

The parameters used in Algorithm [1] are set as follows: 
MIN: minimum cycle length 44 seconds; 
STOPPER: stopper cycle length 64 seconds; 
MAX: maximum cycle length 130 seconds; 
STEP: fixed amount of increment/decrement 6 seconds. 

Fig. IA.141 illustrates the cycle length decision process implemented by Algo- 
rithm [T] The main strategy underlying the above algorithm is to attempt to 
keep the volume ratio within the range [0.85, 0.95]. If the volume ratio was high 
during the previous cycle (over 0.95) the cycle length is increased by a fixed 
amount. Otherwise if green time was wasted on the previous cycle, signaled 
by R < 0.85, the cycle length is decreased. The STOPPER is included to allow 
a steep increase in the cycle length due to a sudden increase in traffic volume, 
when the cycle length is at its minimum. 

The volume ratio of a master node, m, is given by 

R(m) = R(l*,P*), (A.l) 

where I* is the unique inlink flowing in the linked direction, and V* is the linked 
phase. The cycle length of each slave node is equal to that of its master. 

For a non-subsystem node, n, the adaptive cycle length strategy remains 
valid, except that the volume ratio is defined by maximizing R over all inlinks 
and phases, 

R(n) = maxmaxi?(/,P). (A.2) 

Given the cycle length C, the split time of phase V for a master or non- 
subsystem node is given by 

S = r-f^,} . [C — number of phases x S mul — total amber time] + S m - m , (A. 3) 

with the demand function d(V) = max; V(l, V) where I is an inlink of phase V . 
We impose a fixed delay T wait to the nodes each time there is a phase change 
and the next phase does not share any path with the current one. During the 
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Figure A. 14: Cycle length selection by SCATS-like systems. 



time only right-turning vehicles that have been given way to others may traverse 
the intersection. This delay mimics the amber time for phase change. We set 
Tvra.it to 2 second in our simulations, however the precise value does not impact 
greatly on the simulation results provided that the split times are not too small. 
For SCATS, the total amber time in a cycle is 4 second. The minimum split 
time in our simulations was set to S m in — 5 seconds. 

For slave nodes, we demand that the split time of the linked phase must be 
the same as that of its master node. The remaining portion of the cycle is then 
shared between the other phases according to their maximum inlink volumes. 

Initially, at the beginning of each simulation, the cycle length of each node is 
set to the minimum value and the split time plan is chosen based on the turning 
probabilities. We note that since split times are adaptive, the initial condition 
for the splits is unimportant. 

Appendix A. 2. SOTL 

SOTL is an acyclic signal system, in the sense that no fixed ordering of the 
phases is imposed. The following algorithm provides a precise description of 
how SOTL operates at each time step and node. The observable r(n) acts as a 
clock for node n, recording how long the current phase has been activate for. 
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Algorithm 2. SOTL 

Increment r{n) 

for each phase V ^ active ^° 

Increment t(V) 
end for 

if r(n) > S m i n then 

Let it = {v e n = {v u r 2 , ■ ■ ■ } ■ *{v) > 9} 

if n' ^ then 

Lef n" = {Pen': = max P / en / k(V')} 
Let n'" = {Pe n" : t(V) = maxp, £ n" t(P")} 
Uniformly at random choose V G II'" and set V ' ac t{ ve = V 

Set T ( V active) = 
Set r(n) = 
end if 
end if 



When the idle time of node n is greater than the minimal split time, S m i n , 
SOTL chooses the phases for which the threshold functions k are greater than 
the threshold 6, and among those it selects the phases with the maximal k, then 
among those it selects the phases with the longest idle time. If there is more 
than one element in this latter set, then the next active phase will be chosen at 
random from it, however in practice there is typically only one such phase to 
choose from. The fixed delay T wa n is also applied to SOTL. In our simulations, 
the minimal split time was set to ^min = 5 seconds, as was done for SCATS. 
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